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Abstract: We describe new ambitwistor string theories that give rise to the recent 
amplitude formulae for Einstein-Yang-Mills, (Dirac)-Born-Infcld, Galileons and others 
introduced by Cachazo, He and Yuan. In the case of the Einstein-Yang-Mills ampli¬ 
tudes, an important role is played by a novel worldsheet conformal held theory that 
provides the appropriate colour factors precisely without the spurious multitrace terms 
of earlier models that had to be ignored by hand. This is needed to obtain the correct 
multitrace terms that arise when Yang-Mills is coupled to gravity. 
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1 Introduction 


Witten’s twistor-string theory [1] led to remarkably simple formulae for tree-level gauge 
theory amplitudes in four dimensions [2]. Subsequently, the development of new tree- 
level formulae have led the way towards finding the underlying worldsheet models. 
An analogous formula for N = 8 supergravity was discovered [3, 4] and led to the 
construction of a twistor-string theory for N = 8 supergravity [5] . Formulae in arbitrary 
dimensions were then discovered by Cachazo, He and Yuan (CHY) for Einstein and 
Yang-Mills amplitudes [6-8] based on the scattering equations (which underpin all these 
formulae including the original twistor string [9]). These in turn led to the discovery 
of a family of ‘ambitwistor-string’ theories [10, 11]. These are chiral infinite tension 
string theories that provide a natural generalization of twistor-string theories from four 
dimensions to arbitrary dimension. They give the physical theory underlying these 
formulae and lead to new extensions. For example, in the case of the critical type 11 
gravity model in 10 dimensions, the ambitwistor model leads to proposals for how the 
formulae might be extended to loop amplitudes [12-14]. They also give new insights 
into the relationship [15, 16] between asymptotic symmetries and soft theorems [17-20]. 

Recently, CHY discovered new formulae for a large collection of theories. These 
include the original Einstein (E), Yang-Mills (YM) and biadjoint scalar (BS), together 
with new formulae for Einstein Maxwell (EM), Einstein Yang-Mills (EYM), (Dirac)- 
Born-Infeld ((D)BI), Galileons (G), Yang-Mills Scalar (YMS) and nonlinear sigma 
model (NLSM) [21, 22], see figure 1 below for their diagram 1 of new amplitude formu¬ 
lae and relationships between them. They raised the challenge to find the underlying 
ambitwistor string theories that give rise to these formulae. Indeed, Ohmori in parallel 
work has already found the ambitwistor strings for the BI and Galileon theories [23]. 
It is the purpose of this paper to review and to give further details for his results and 
to show how the remaining new CHY formulae also arise from new ambitwistor string 
models. 

The formulae of Cachazo, He and Yuan give rise to formulae for scattering am¬ 
plitudes as a sum over solutions to the scattering equations of certain Pfaffian and 
determinant expressions. They give tree amplitudes in the general form of an integral 
over n copies of the Riemann sphere 



where a t , i = 1 ,,n are complex coordinates on each of the Riemann spheres, ki the 

1 We thank CHY for permission to reproduce their diagram. 
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Figure 1. Theories studied by CHY and operations relating them. 


null momenta of the massless particles in the scattering process, and 


P“(a) = Y, 

i=l 


K 

cr - (Ti 
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and G = SX(2, C) x C 3 is the residual gauge symmetry of the ambitwistor string 
fixed according to the standard Fadeev-Popov procedure. The integrand naturally 
decomposes into factors I 1 and I r that depend on the a,;, ki and the polarization and/or 
colour data of the particles whose scattering is being computed and depends on the 
theory. The delta functions 

S(z) = 3- —— = <J(9te)<5(3f2)d;z, 

Z7 TIZ 

where 3ft denotes the real part and 9 the imaginary part, impose the scattering equa¬ 
tions: ki ■ P(cri) = 0. The integrals essentially reduce to a sum over solutions to the 
scattering equations of I l I r multiplied by a Jacobian factor. The I 1 and I r can be cho¬ 
sen from five different choices and the various theories arise from the different possible 
combinations. 
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Table 1 . Theories arising from the different choices of matter models. 


Ambitwistor strings are chiral infinite tension analogues of RNS strings that can 
be interpreted, after reduction of constraints, as strings whose target space is the space 
of complexified null geodesics in Minkowski space. This space of complexified null 
geodesics has become known as ambitwistor space. Ambitwistor strings are built out 
of a basic bosonic model together with worldsheet matter. The bosonic model leads 
to a framework in which the vertex operators required for amplitude calculations in¬ 
corporate the scattering equations. The vertex operators also allow for the insertion of 
two currents v l and v r and these can be constructed from additional worldsheet matter 
(the natural choice for v l and v r in the bosonic model does not seem to lead to inter¬ 
esting amplitudes). The various Pfaffians, determinants or Parke-Taylor factors that 
are possible choices for the I 1 and P arise as worldsheet correlators of currents for the 
v l and v r respectively. Corresponding to the five choices for the I 1 and P in the CHY 
formulae we will introduce five choices of worldsheet matter, see table 1. 

In the original models of [10] just two ingredients were used to construct I 1 and 
J r , worldsheet supersymmetry ST, and a current algebra Sj. Einstein, Yang-Mills and 
Biadjoint scalar theories were obtained from the choices ( S l ,S r ) = (ST, ST), (ST, ST) 
and (Sj, Sj ) respectively. The current algebra Sj has the defect that it also leads to 
multi-trace terms in its correlators that were ignored by hand. Here we use a different 
worldsheet CFT, the comb system 2 , Scs- This gives a new way to obtain colour factors 

2 This was originally introduced by David Skinner and one of us [24] in the context of twistor-strings, 
but never published. 
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together with their Parke-Taylor cyclic denominators in such a way that these multi¬ 
trace terms simply do not appear. Furthermore, the colour factors are presented not 
as cyclic single trace terms, but as strings of structure constants arranged in a ‘comb’, 
hence the name. However, the number of gauge particles in this system is doubled. To 
remedy this issue, a reduced system Sym with the correct number of gauge particles can 
be constructed, but this is is always anomalous. Nevertheless, it is sufficient to produce 
the correct tree amplitudes and so we use this system instead of the current algebra in 
the table 1. It can be replaced by Scs if we are seeking an anomaly-free theory, but 
then we must accept the doubling of gauge particles. The remaining systems that we 
use will be combinations of these (with S p ^ essentially being the abelian limit of the 
combination Scs.'V of the comb system with worldsheet supersymmetry). These will 
be described in more detail in the main text. 

There are a number of questions that one can ask about these models. For example, 
if they are critical and anomaly free, then one can attempt to calculate loop amplitudes 
by taking the correlation functions on higher genus Riemann surfaces as described in 
[12]. For this to work at 1-loop, we must check modularity. Another issue is as to 
whether there are any further vertex operators in the theories and if so, we can hope 
to extend the theory to include additional fields and calculate the corresponding am¬ 
plitudes. This in particular happens in theories containing worldsheet supersymmetry 
Sy and leads to supersymmetric extensions of the theories and amplitudes as described 
also in [12]. We will find a number of new critical theories and give a brief discussion 
of these issues in the conclusions section. 

Potentially the most interesting of these models is that for Einstein-Yang-Mills. 
We obtain these in two forms. One gives the correct CHY tree-level amplitudes, but is 
anomalous using Sym- The other has vanishing central charge in 10 dimensions, but 
has doubled gluons in the theory. The gauge theory part of the action is given by 



( 1 . 1 ) 


and we refer to it as T*YM as it describes a linearized Yang-Mills field a propagating 
on a full Yang-Mills background for the field A with curvature F. Here a is canonically 


conjugate to F hence the name T*YM as opposed to TYM. This should give correct 


Yang-Mills amplitudes at one loop but has no higher loop amplitudes in the pure gauge 
sector. In its critical dimension d = 10, we would expect it to give a valid expression 
for the 1-loop integrand for Yang-Mills also. 

Table 1, showing how the theories are determined in terms of a pair of worldsheet 
systems, is a remarkable manifestation of the notion of double copy. This notion has 
been explored mostly in the context of gravity amplitudes, which are obtained as the 
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double copy of gauge theory ones [25, 26]. In the formalism of the scattering equations, 
this is the double copy of Pfaffian factors, and in ambitwistor string theory, this is the 
double copy of the worldsheet system Sy, as in table 1. The amplitude formulae of 
ref. [22] and our results extend this notion to a range of other theories. Regarding 
the relation to previous work, we should mention that a double copy construction for 
Einstein-Yang-Mills amplitudes was first presented in [27] for the single trace contri¬ 
bution, and in [28] for the complete amplitude, with results also at loop level. These 
double copy constructions are based on the colour-kinematics duality [26, 29], whose 
relation to the scattering equations has been explored in [6, 30, 31]. 

2 Ambitwistor string models 

2.1 The bosonic ambitwistor string 

The original bosonic ambitwistor string starts with an action obtained by complexifying 
the standard massless superparticle. We complexify the worldline to be a Riemann 
surface £ with holomorphic coordinate a G C, and the target space is taken to be the 
cotangent bundle T*C d of complexified space-time (so that the (P /i; X M ), n, v — 1,... d 
are holomorphic coordinates and space-time has a holomorphic flat metric r/ /w ). We 
use the bosonic action 

S B = S B [X, P] = i- J P ■ BX + eP • P, (2.1) 

where dX = dd d&X and the will be used to denote contraction of indices whenever 
the contraction is unambiguous. We take the components of P p to be complex (1,0)- 
forms on the worldsheet, so that suppressing space-time indices, P = P CT (a) da G 
K := T*£. It follows that e must be a (0,l)-form on £ with values in T := T£ 
the holomorphic tangent bundle. It plays the role both of a Lagrange multiplier that 
enforces P 2 = 0 and of a gauge field for the transformation 

be = dot , 5X = aP, 5P = 0 . 

Restricting to P 2 = 0 and reducing by the gauge freedom reduces us to ambitwistor 
space, the space of complex null geodesics. 

This gauge freedom can be fixed by setting e = 0 with the introduction of ghosts 
(b,c) G (Ji 2 ,T), where K = T*£ is the bundle of (l,0)-forms. The usual diffeomor- 
phism freedom can also be parametrized using d e = 5 + ed and fixed with e = 0, 
together with the introduction of the usual ghosts (6, c) G (Jl 2 ,T). This leads to the 
basic BRST operator 

Qb — / cT + cP 2 . 


6 


It was shown in [10] (see also [12] for a more systematic treatment) that these 
structures are sufficient to lead to a framework in which amplitudes are computed 
using vertex operators built from some operator V G K 2 . This will be the case for all 
our ambitwistor-string theories. When E = CP 1 , there are three zero-modes each for c 
and c, and three vertex operators are fixed, coming as CjCjV)e* fei ' x at cq, i — 1,2,3. The 
remaining n — 3 are integrated 

V* = f • P(ai))e iki ' x Vi(ai ), 

where 

P{&) = V —, 8{z) = = S(^z)6(^z)dz. 

U — Uj 27 TIZ 

j J 

The delta-functions 5(ki- P(cq)) impose the scattering equations that provide the back¬ 
bone of all the CHY formulae and their twistor-string precursors. We will take these 
aspects of the scattering amplitude calculations for granted in the following and will 
not mention them further. 

2.2 Vertex operators and worldsheet matter 

In general, we will consider theories with actions of the form 

s B + s l + s r 

where S l and S r are distinct matter theories on E that will be described in more detail 
in the next two sections. These will contribute to our formulae for the vertex operators 
V G K 2 . The vertex operators all have an e lk ' x factor with the remainder V factorizing 
into two independent currents, 

V = v l v\ v\v r eK. 

The v l , v r will be constructed from the matter models, S l and S r respectively and 
constrained by quantum consistency, BRST invariance, and perhaps further discrete 
symmetries. Invariance under Q B for example implies k 2 = 0 because the P 2 term in 
the BRST operator brings down k 2 in its double contraction with e lk ' x . We will also 
use the notation u l , u r for such currents when they are fixed with respect to fermionic 
symmetries, see below. 

Essentially the only candidate for v l and v r in the purely bosonic model above is 
e-P for some polarization vector defined up to multiples of k M under Q B equivalence. 
This leads to unphysical formulae for gravity amplitudes, or at least with no clear in¬ 
terpretation. In order to obtain more interesting models, we will introduce worldsheet 
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matter models S l and S r that will generate the currents v l and v r in the vertex oper¬ 
ators. In general, we will take the models S l and S r to be distinct matter theories so 
that the correlator will factorize into a product of one for the left currents and one for 
the right currents, and we will be able to calculate them separately. In order to ensure 
that the only allowed vertex operators do indeed factorize in this way, we will impose 
discrete symmetries that are analogues of the GSO symmetries of conventional string 
theories, and we will use that name for these symmetries as well. 

3 Worldsheet matter models and their correlators 

In [10], two matter models were considered: (1) S p , a current algebra which we will take 
to be generated by free fermions, and (2) Sy, which introduces a degenerate worldsheet 
supersymmetry. This latter extends Q so as to change the choice of current v = e ■ P 
in the bosonic model to one that we will want. These led to three models with (S l , S r ) 
given by (S x y l , SV r .) for type II supergravity, (S^,S p ) for Yang-Mills amplitudes and 
(S pi ,S Pr ) for amplitudes of a biadjoint scalar theory. In this paper we will consider a 
third type of matter that we call the ‘comb system’ Scs, [24], a worldsheet conformal 
held theory that will be important for Yang-Mills amplitudes so called because its 
correlators give colour invariants in the form of comb structures built out of structure 
constants rather than colour traces. In the rest of this section, we describe these matter 
systems, and the natural currents to which they give rise as candidates for v l and v r 
and their correlation functions. In the next section we see how these are altered when 
these systems are combined. 

3.1 Free fermions S p and current algebras Sj. 

The standard action for ‘real’ free fermions p a G A' 1 / 2 , a — 1,... m, is 

S„= f p‘Bp° , 

(the summation convention is assumed). The term ‘real’ is used to distinguish them 
from the complex fermion system given by 

J ~ paBpa ' 

The simplest currents in the real case are j ab = p a p b and form an elementary example of 
a current algebra for SO (in ) (in the complex case = pbp a generate a current algebra 
for SU(m)). 


More generally, we can consider an arbitrary current algebra j a G K Cg> 0 , where 0 
is some Lie algebra, a — 1,, dini 0 , satisfying the usual current algebra OPE, 


j»A 0) 




ab 




+ 


if abc t c (a) 


a 


+ ... , 


(3.1) 


where f abc are the structure coefficients, [t a ,t h ] = f abc t c , S ab is the Killing form, and k 
is the level. This could be contracted from free fermions, WZW models or some other 
construction and we will generally represent such matter as S 3 . 

Given choices of t G 0 , the current algebra can contribute 


v — t-j 


to one or both factors Vi and v r of the vertex operators V. The current correlators 
( t\ ■ j i... t n ■ j n ), where j t = j(< 7 j), lead to Park-Taylor factors: 


PT(l,...,n) 


trgi ...t n ) 

°'120‘23 • • • °nl 


where — a r However, the correlators also lead to multi-trace terms that are 

ultimately problematic and unwanted. 


3.2 Worldsheet suspersymmetry S^. 

Worldsheet supersymmetry is introduced by adding fermionic worldsheet spinor holds 
vE' g C d ®K 1 / 2 , and a gauge held x G n(°’ 1 )(T 1,/2 ) for the supersymmetry. Their action 
is 

The constraint leads to worldsheet gauge transformations 


j T • <9T + yP • T . 


5 X = dr ], SX = /;M/, hT = r]P , 8P = 0 , 

where 7] is a fermionic parameter. Gauge fixing leads to bosonic ghosts 7 G T 1 ^' 2 and 
corresponding antighosts f3. The BRST operator acquires an extra term 

Q\v = jl 'yG'S ’, Gy := P • T . 

On CP 1 , the ghosts 7 have two zero modes. Thus, as far as the fermionic symmetry is 
concerned, we need two hxed vertex operators with one current factor of the form 5 ( 7 ) 
multiplied by a held now with values in A' 1 / 2 , and then the ‘integrated’ ones (in the 
fermionic sense) arising from descent. The relevant currents are 


11 = 5(7)e - T, u = e- P + /c-Te-T, 
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with just two of the us required in a correlator. 

These operators are invariant under the discrete symmetry that changes the sign 
of T, x and the ghosts. Imposing invariance under this symmetry will exclude mixing 
between the ingredients of these operators thought of as parts of S l and others that 
might be part of S r . We will refer to this as GSO symmetry. 

The correlators of these currents lead to the reduced Pfafhans of CHY: 

(■ uiu 2 v 3 ... v n ) = Pf'(M) = ---Pf(Mi 2 ), 

cti — cr 2 

where M is the skew 2 n x 2 n matrix with n x n block decomposition 


( \ 

A -C T 


Aij 


ki ■ kj 


a 


y 


Bij 


d • e j 
°ij 


and 

Qj = , i 7 ^ j, Ca = -Ci ■ P{ai ), 

and M 12 is M with the hrst two rows and columns removed. 


3.3 Comb system Scs- 

The comb system [24] was introduced as a way of obtaining colour factors as sequences 
of contractions of structure constants rather than as colour ordered traces. In general, 
such contractions can be generated from trivalent diagrams with the structure constants 
f abc of some Lie algebra at the vertices and contractions 5 ab along the internal edges. 
It is well known that these are linearly dependent as a consequence of the Klciss-Kuijf 
relations with a basis being given by ‘combs’, with n — 2 vertices lined up in a row 
[32, 33] and end points given by 1 and n: 


2 3 

1 —I—L 




j^a 1 a 2 b 1 J^b 1 a 3 b 2 . . . y6 n _3a n _ia„ 


The comb system has the remarkable property that, in conjunction with worldsheet 
supersymmetry, only these combs arise from correlators and not the multitrace terms 
that arise from an ordinary current algebra. This system arises from an action for 
matter fields p,p,q,y G g ® A ' 1 / 2 i.e., worldsheet spinors taking values in the Lie 
algebra g of some gauge group. The worldsheet action is 

Scs = / P ■ dp + q ■ dy + XP ■ Q[p, P] + [q, 2 /]) 
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with p, p fermionic and q, y bosonic and the • is used to denote the Killing form on the 
Lie algebra. As before, x is a gauge held on the worldsheet with values in T 1 / 2 <g) fi 0,1 
and we are gauging the current 3 p ■ (^[p, p] + [ q,y ]), which is a section of K 3 / 2 . The 
gauging introduces transformations now for fermionic a G T 1//2 

S{p,p,q,y) = a (\lP,Pl \P,P\, \p,q\, [PiV]) , <^X = <9a. 

As in the case of worldsheet supersymmetry, gauge fixing gives bosonic ghosts 7 G T 1 / 2 
and antighosts (3 with a contribution to the BRST operator of 

Qcs = j> iQcs , Gcs ■= P ■ Q[p, p] + [g, ?/]) • 

As for Sq,, there are two zero-modes for the ghosts, and so we will need two fermioni- 
cally hxed operators with the rest integrated. The currents that contribute to the vertex 
operators in this system now depend on a Lie algebra element t G 0, with two types of 
hxed and integrated ones respectively being 

u = S{'y)t-p, u = 8 {i)t-p, v=-t-\p,p], v = t- ([p,p] + [q,l]) . 


Here v = {Qcs,u} and v = {Qcs,u}, and, in any correlator, we need two hxed and 
the remaining unfixed vertex operators 4 . Notice that v — t ■ j, where j a is a level zero 
current algebra, and that 


v(cr) t' ■ p( 0 ) ~ 


it, ?} • m + 

cr 


v(a) t' ■ p( 0 ) ~ 


[t, t'} ■ P(0) + 
a 


The correlators are as follows 


Proposition 3.1 (Casali-Skinner) Correlators of the currents u, v, u, v are only non¬ 
vanishing when there is just one untilded current and give 

{U1V2 ■ ■ ■ Vn-lUn) = (U1V2V3 ■ ■ ■ V n -1 u n ) = C(l, . . . , 7 l) 

3 With different assignment of worldsheet spins this current would be a normal BRST current. If 
we were to take p, y , scalars and p, q sections of K , then p and p could be taken to be the ghosts 
associated to gauge fixing a worldsheet gauge field a £ fl 0,1 1S> g with action f s q ■ By + q ■ [a,y). This 
fact allows one to see the consistency of this current reasonably rapidly. 

4 A more symmetric way to understand this is to say that we choose all unintegrated vertex 
operators, but then we must insert n — 2 ‘pi c t ure "Changing operators’ 

T = 5(/3)p- Q[p,p] + [g,y]) • 

These could be inserted anywhere in general. If inserted at one of the u, u insertion points, it will 
convert it into a corresponding v,v. A similar approach can be taken for correlators associated with 
the Sq, matter system. 
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where 


C n = C{ 1 , ...,n) 


tr(ti[t 2 , [..., [t ra -i,t ra ] ■ ■ ■]]) 

012023 • • • &nl 


+ Perm(2,..., n — 1 ). 


Instead of giving the colour traces, we obtain ‘combs’, i.e., strings of structure 
constants tr(fi[t 2 , [..., [t n _i,t n ]...]]) as described in [33, 34], 

The argument is as follows. The fact that we can have at most two u, us is the 
standard counting of 7 ghost zero modes. Consider the p, p contractions: that these 
are the only nontrivial correlators comes from the need to have as many ps as p, so it is 
easily seen that we can have only one untilded current which can either be a u or a v. 
The v,vs connect along a ‘comb’, whereas the u,us form the ends. Such contractions 
connecting all n vertex operators form the right hand side above. We can also have 
contractions in which a collection of vs come together in contractions to form a loop. 
This is where the (q, y) system comes into play. These can only form loops, but, being 
bosonic, their loop contractions cancel such loop contractions from the p,p system. 
This can also be seen from the form of the current algebra generated by the vs. This 
has by construction level zero so that, after a sequence of OPE’s, cannot generate a 
nontrivial trace. 


3.4 Other systems with comb structure, Sy m■ 

A problem with the CS system above is that there are clearly two types of gluon, tilded 
und untilded corresponding to the vertex operators (u, v) and (u, v) respectively. We 
will see that this is not appropriate for pure Yang-Mills although it does give a theory 
that is sufficient to generate Einstein-YM tree amplitudes correctly on certain trace 
sectors, the ones selected by the choice of untilded operators . 5 The system we introduce 
here will give the complete Einstein Yang-Mills amplitude from a single correlator, but 
will be anomalous. 

A worldsheet CFT that will generate YM following the ideas above requires the 
following ingredients. We need a fermionic worldsheet spinor p a E g for the fixed vertex 
operator, a current algebra v a E 0 at level zero for the integrated one; the level zero 
allows us to avoid multitrace terms and loops. Finally we need a spin 3/2 current Gym 
with the following OPE to give the appropriate group compatibilities and descent: 


raft 

P»/(0) ~ — , v a (a)p\ 0) 
a 


f abc p c ( 0) 


cr 


£(a)/P(0) 


u“(0) 


a 


G(&)G ( 0 ) 0 . 

(3,2) 


5 One may try to symmetrise the correlator in tilded versus untilded gluonic operators, for instance 
by using u t + Ut and v t + v t . but then there will be an over-counting of contributions, so that the 
relative factors of different terms are not correct. 
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It is easy to see that this can be partially realized with p a a ‘real’ free fermion with 
action | f p a dp a and with 

va = ~\f abc P b P c + j a , j»/( o) ~ o, 

we will obtain the first two of the equations above. In order for v a to be a current 
algebra with level zero, because \f abc p b p c is a current algebra with level — C where 
j'abc fobc _ we mus t t a ke ja p e a curren t algebra with level k = C. There are 

many ways to do this, so let us leave this to one side for a moment. We then need to 
construct Q. In order for Q to generate v a from p a , we must have 

G = -\f abc p a p b p c + p a j a + ■■■ 

D 

where the ... has nonsingular OPE with p a and j a . At this point, however, we see that 
an anomaly arises preventing {G,G} = 0. To be specific, 


. C'dim(G) ■ j a j a ( 0) : 

G(a)G 0- sU. + J l AT 

a A a 

where we recall that the energy-momentum tensor of the current algebra j is given by 
T((j) =: j a j a (cr) ■ /2k. Therefore, we are able to satisfy the first three equations of 
(3.2), while the last equation is anomalous. 


3.5 Central charges 

We remark that the theories Sb, S p , Sy and Scs above respectively have central charges 

cb = 2d — 52, c p = m/2, = d/ 2 + 11, ccs — 11, 

the latter being just that of the (3 — 7 system as the dim G parts cancel via supersym¬ 
metry. (This can be different if the (q,y) are not taken to be spin 1/2.) Notably, the 
type II supergravity model is critical in 10 dimensions as then eg + 2c^, = 0. These 
considerations are less interesting for Sym as that theory is already quite anomalous, 
and in any case its central charge will depend on the choice of current algebra j a . 


4 Combined Matter models 

On their own, the new worldsheet matter theories Scs and Sym of the previous section 
do little more than give an alternative to the current algebras in the original models 
of [10] that avoids the multitrace terms that were neglected by hand. To obtain new 
theories, we will consider the contributions to S l or S r of combinations of the above 
matter systems. Even without Scs and Sym, we will obtain a number of interesting 
new models. Here we will consider the allowable vertex operators and the correlators 
of the various combinations that we can form. These are summarized in the table 2. 
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Fermionic current Q 

Matter 

Vertex operators 

Correlator 

sv 

P 4 - 


uy = 5 ( 7) e • \I/ 

Pf(M) 


P 

P-* 2 

^1, <f 2 

uq, 1 = 5(72) k ■ 'Fi 

uq, 2 = < 5 ( 7 i) k ■ 'F 2 

(Pf(A)) 2 


P 4 - 

p a a = 1,... ,m 

= (ify) e • 

U Pa = Kl) Pa 

Pf( X )Pf(M| red ) 

Scs 

P ■ '&+tr(p(±[p,p\ + [q.y])) 

^ (A p), {q,y) 

u^t = 5 ( 7) e • \I/ 

ucs = 8(7) tr (tp) 

ucs = £(7) tr {tp) 

C (1) ...C (m) Pf(n) 

Scs 

t r {p(\[p,p\ + [q,y])) 

( p , p ), (q,y) 

ucs = £(7) tr (tp) 

ucs = £(7) tr {tp) 

C n 


Table 2. Table of matter models, their combinations and worldsheet correlators 


4.1 S P) qj 
Here we take 

S l = S Pt * := S p + S*. 

Although the free fermion system S p would seem to naturally lead to the SO(m ) current 
algebra j ab = j ^ = p a p b 1 and therefore superficially be thought to give the same results 
as the current algebra, in the presence of worldsheet supersymmetry, the currents j ab 
as constituents of vertex operators are not BRST invariant, since 

{<2*, j ab e ik X } = ik ■ \I>j ab e tk ' x ± 0 . 

However, in this context allowable fixed and integrated currents are respectively 

u a = 5(~f)p\ v a = k-^p a , a — 1,..., m. 

We also have the standard BRST invariant currents from S which in this context we 
will denote u € = <5(y)e ■ T and v e = e ■ P + k ■ Te • T. 

In general we will be concerned with a correlator ( 111 U 2 V 3 ... v n ) where, if (7, h) is a 
partition of 1 ,..., n, for * £ 7 the current will be one of the new photon currents, and 
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for i E h it will be a Sq, current depending on a polarization vector e M . The correlator 
will factorize into one for the constituent p a s and one for the Ts. We compute these 
as Pfaffians of the associated matrices of possible contractions in the correlator. The 
simplest is the p system. If we restrict it to take values in an algebra with vanishing 
structure constants, e.g. © m u(l), the OPEs lead to the |y| x |y| CHY matrix 

$ a i a i 

Xij = -, i,jE 7, i ^ j, otherwise X l3 = 0. 

a ij 

The Kronecker delta tr (t a H aj ) in the numerator ensures only photons of the same flavour 
interact. 

Much as before, the T system leads to the matrix of possible T contractions 


^ A 11 A lh 


MRed ~ 


7 J^hh _ ^jhhT 


C hh B 


where we have divided the matrix into the bock decomposition under n — |^y| + \ h\ and 
A^ = — , i -f- ji An = 0 , Bij = - - , i, j E h,i 7^ 3- 

(Jij (Jij 

and 

' kj ■ ■ . 

Cij — , i E h,i fz j . 

Finally, the additional e • P term in the Sq, vertex operator is incorporated by setting 
Cu = —€i ■ P(<Ji) as before. In this case, we obtain a reduced Pfaffian associated with 
the two fixed vertex operators as before. Our final correlator expression is therefore 


( Ul u 2 v 3 ...v n )= Pf(Y)Pf \M Red ). 


Now for the GSO symmetry we require all fields, p, T and the ghosts to change 
sign simultaneously. 


4.2 Sq,^yq, 2 

Here we take two worldsheet supersymmetries 

s l = Sq,, + Sq, 2 . 
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There are two contributions to the BRST operator, Q^ 1 + Qq, 2 . The normal currents 
from and Sy 2 are no longer invariant as, for example, 


Q* 2 < 5 ( 7l )e • T ie ife ' x = 5 ( 7l )fc ■ T 2 e • T ie ifc ' x + 0. 


However, the nontrivial BRST invariant currents are descendants simply of £( 71 )<5( 72 ) 
so 

u = 5 ( 71 ) 5 ( 72 ), v = k • dRfc • dR, 

as given in [23] (and we also have partial descendants 5 ( 7 i)fc- d / 2 and ( 5 ( 7 2 )fc-'I / 1 ). Again, 
the correlator of n such vertex operators factorizes into a product of the Pfaffians of the 
matrix of all possible dR contractions and that for all dR contractions. These matrices 
are given simply by the A matrix with off-diagonal entries ki ■ kj/ctij as before. This 
matrix has co-rank two and we take a reduced Pfaffian (corresponding to the choice of 
fixed versus integrated vertex operators). We therefore now obtain 


(U 1 U 2 U 3 U 4 • • • v n ) = Pf'(A ) 2 . 


One might ask whether one can carry on to combine three or more SR systems, but 
this is not possible on one side, that is, to produce nontrivial BRST invariant currents. 

Again for the GSO symmetry we require all fields, dR, dR and the ghosts to change 
sign simultaneously. 


4.3 Sym,v 

In Sections 3.3 and 3.4, we introduced Scs an d Sym whose correlators provide the 


colour comb-structure together with Parke-Taylor factors. For the remainder of this 
section, we will combine each of these two systems with S 4 ,. The goal is to obtain the 
building block of Einstein-Yang-Mills amplitudes that gives the appropriate interactions 
between gluons and gravitons. We start by discussing the combined theory Sym,v, 
which is slightly simpler than Scs and possesses the main important features. Despite 
Sym not being quantum-mechanically consistent - and this problem extends to Sym,v 
- we are able to obtain tree amplitudes. The theory Scs,v can be made consistent, but 
has two types of gluons and the corresponding amplitudes arise from an action that is 
not Yang-Mills (although it contains its classical solutions). 

Since both worldsheet matter theories Sy and Sym involve the gauging of spin 3/2 
currents Qq, = P ■ T and Gym = P • (—|[p, p\ + j ), we have the option of gauging both 
these currents together or separately. If we gauge them separately, we find that the 
resulting system is too restrictive to lead to interesting results. Thus we gauge the sum 
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perform gauge fixing and introduce a single set of ghosts (/3, 7 ). We find that the 
currents 

u t — 8{l)p ■ t , u e = 5{ 

still give us allowed hxed vertex operators. BRST descent leads to the integrated vertex 
operators 

v t = k • Tp ■ t + , v e — e ■ P + e ■ ■ 4/, 

where denotes the original Sym integrated vertex operator, satisfying the OPE 
relations (3.2) except the last. Although the failure of the last relation means that 
the BRST quantisation is inconsistent, the correlator of the vertex operators does 
nevertheless give the correct amplitudes. 

In the previous section, we saw that the system Sym on its own gives the correct 
colour-dressed Parke-Taylor factors, in terms of a comb structure. The combination 
with Sq, leads to additional insertions of p ■ t and these will start additional combs. In 
this way we obtain multiple colour combs/traces and get the right interactions with 
gravity states. On the other hand, the system Sty on its own leads to a reduced 
Pfafhan. The combination with Sym will lead to a different but closely related Pfaffian 
that incorporates the multi-comb structure. We now describe the complete correlator. 

Theorem 1 As in [22], let the sets g index the gluons with vertex operators u t ,v t , and 
h the gravitons with vertex operators u e ,v e . To be non-zero, a correlator must contain 
two fixed vertex operators u’s, with the remaining ones being v’s. The correlator is then 
a sum over all partitions of the gluons into sets Ti,T 2 ,... ,T m , where U^ = 1 T i = g and 
|Tj| > 2. Each partition gives rise to the term 



( 

\ 


A-ab ^4 acj -A a dj 

( -C T ) ab 


■4 ab A-acj A-Cidj 

(~c T U 

/C( Cl , di\Tf) ■ ■ ■ JC(c n , d n \T n ) Pf 



ci<diSTi 

A-dib -4 diCj A-didj 

(- c T \ kb 





Cab C acj C adj 

B a b y 


Here, a,b label gravitons and Ci,di label gluons in T t , so that A a b is an \h\ x \h\ matrix, 
A c . b is an m x \h\ matrix, and A c . c . is an m x m matrix. Moreover, we defined 

lC(i,j\T) = <r ji C(T), (4.2) 

where C(T) is C n restricted to g e T. The reduced Pfaffian Pf is defined in eq. (A.15). 
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The proof is given in appendix A. This correlator reproduces the main building 
block of the CHY formula for Einstein-Yang-Mills amplitudes in [22], Although not 
quite in the same form, the equivalence can easily be seen from Eqs. (3.16) and (3.17) 
of [ 22 ] and this form is more natural from its derivation as a correlator. 

4.4 S C s,<s> 

While Sym,v gives the correct amplitude, its BRST quantisation is inconsistent. We 
can obtain the same structure from Scs ,9 by combining the worldsheet theories Sy and 
Scs 5 which has the advantage of being anomaly free but the disadvantage of containing 
two types of gluons. 

As for Sym.'V we gauge the sum of spin 3/2 currents = P ■ T and Qcs = 
p ■ (f [p, p] + [q,y]), introducing the action 

Scs* = f 'h • dv + p ■ Bp + q ■ By + x (p ■ ^ + P ■ (^[p,p] + [g, 2/])) • 

Now the Lie-algebra valued fermion p is complex (i.e., not equal to p), unlike the 
previous case of Sym,v- This will change the physical content of the model. The gauge 
fixing of y introduces just one set of ghosts (^, 7 ), and we find the standard fixed 
currents for Scs an d Sy, 

u t — S(pf)p ■ 1 1 u t — 5 (p/)p ■ t , u e = <5(7)e •'L. 

The BRST descent leads to the following currents 

v t — k ■ • t + , v t — k ■ typ ■ t + , V £ — e ■ P + e ■ ^>k ■ \k, 

where u/ and 7/ denote the original Scs integrated vertex operators, so that v t and v t 
acquire a new term in T. 

To impose GSO symmetry, we require invariance under flipping the sign of the 
fields p, p, q, y, \H, x an d the corresponding ghosts. 

Since we have untilded vertex operators u t ,v t , and tilded ones u t) v t) the correlator 
will depend not only on the number of gluonic vertex operators versus gravity ones 
u e ,v e , but also on the choice of whether the gluonic operators are of untilded or tilded 
type. Recall from the previous section that, for the theory Scs on its own, the only 
non-vanishing correlators were those with a single untilded operator and this led to a 
single comb colour structure that is equivalent to a single trace term. This followed 
because of the need to have the same number of ps and ps in a nontrivial correlator 
and a single p could only arise in one or both of the two fixed vertex operator. Now 
single ps appear in v t and this essentially represents the coupling to gravity. Thus the 
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coupling to gravity introduces multiple trace terms, with the interaction between each 
single trace structure being mediated by gravity. It is easy to see that with the Scs 
system we can now have as many untilded vertex operators as we like with their number 
corresponding precisely to the number of traces. 

Theorem 2 Let the set g index the gluons and h the gravitons. To be non-vanishing , 
a Scs,v correlator must have two fixed vertex operators, with the remaining ones inte¬ 
grated. The correlator of such a collection of vertex operators is a sum over all partitions 
of the gluons into sets Tf, T 2 ,..., T rn , where m is the number of untilded gluonic vertex 
operators, and such that there is only one such vertex operator per T i} U™ jTj = g, 
\Tf\ > 2. Each allowed partition gives a contribution equal to (4.1). 

Thus the correlator is the same as for Sym,v, except that there is a restriction on 
the allowed partitions of the gluons into traces. 

5 New Ambitwistor String Theories 

We can now assemble the full table of theories by combining the various possible choices 
of matter models on the left and right. These can be identified with their corresponding 
space-time theories by comparing the correlators to the formulae of CHY, and this 
results in table 3. Hopefully the acronyms for the models are self-explanatory except 
perhaps that BS denotes the bi-adjoint scalar <f aa ' , where a and a' are respectively 
indices for the Lie algebras of SU(iV) and SU(iV 7 ), with action 



where f abc and f a ' b ' c ' are the structure constants of SU(iV) and SU(AT / ) respectively. 
Galileon theories are described by the action 



where g m are freely prescribable parameters. However, our amplitudes only have one 
parameter. The theory that is singled out by our model is the one described in [35] in 
four dimensions, and in [22, 36] in general dimension, which has smoother soft behaviour 
than the generic ones. 

The Born-Infeld action is 
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\ s r 
s l \ 

S* 

Sh-i ,^2 

c(m') 

q(N') 

°YM,V 

o(N') 

°YM 

s * 

E 





Sy 1,^2 

BI 

Galileon 




qM 

h^lu(i) m 

DBI 

EMS| U ( 1 ) m0U ( 1 ^ m 



o(N) 

°YM 

EYM 

extended DBI 

EYMS | SU ( A T) (g)U ( 1 )m' 

EYMS | SU ( iV ) ( g )SU ( A r/) 


oW 

°YM 

YM 

NLSM 

YMS|g U ( JV ) 0 u(i)m' 

gen. YMS | (jv') 

BS 


Table 3. Theories arising from the different choices of matter models. 


the Dirac Born-Infeld is 

Sdbi = J d D x ^ ^y/-det - tF - 1^ , 

and the nonlinear-sigma model is 

SWlsm = J d D x f-^tr ((1 - A 2 $)- 1 ^$(l - A 2 <3>) _1 <9 M <I>) ^ , where $ = . 

In the table 3, we have only used Sym■ Although this is sufficient to produce the 
correct tree-level amplitudes, it is an anomalous matter system and so has no hope 
to be extended beyond tree-level, and indeed its meaning as a string theory is unclear 
even at tree level. We can obtain the same tree-amplitudes up to combinatorial factors 
by use of the comb system Scs and this is not anomalous. However, this does lead to a 
doubling of the gauge degrees of freedom as described below in detail for the Einstein 
Yang-Mills system and bi-adjoint scalar. 

In table 4 we list the vertex operators in each model and the central charges. Setting 
the central charge to zero gives the models that are critical and for which there is some 
reasonable hope that loop integrands can be described via these theories if they prove 
to be modular. 

5.1 Einstein Yang-Mills and T*YM 

The worldsheet model that we discussed in the context of Einstein Yang-Mills theory, 
Scs,^ i has a consistent quantisation. On the other hand, it does not correspond strictly 
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to the building block of Einstein-Yang-Mills amplitudes, because only trace/comb struc¬ 
tures consistent with the choice of untilded vertex operators are allowed. Attempts to 
End a theory that reproduces this correlator seem to lead back to the anomalous Sym,v 
system. 

Since the theory Scs,v presents no problems, and has correlators which match part 
of the Einstein-Yang-Mills building block, it is natural to ask whether it is related to 
a known theory. This theory must contain two types of gluons, associated to tilded 
and untilded vertex operators, and the untilded type must give the number of allowed 
multiple trace terms in an amplitude. These conditions are satisfied by the following 
spacetime action for the gauge field 



(5.1) 


The field a^ is a Lagrange multiplier enforcing the Yang-Mills equations, D U F^ U = 0, 
and the action can be seen as a linearisation of the Yang-Mills action, + a M . 

The field A M corresponds to the tilded degrees of freedom, and the field corresponds 
to the untilded ones. Since the propagator of this action connects to A M and the 
vertices contain a single a M , the Feynman rules and a straightforward graph-theoretic 
argument show that there is one and only one external field per trace, also when 
the system is minimally coupled to gravity. 

5.2 Bi-adjoint scalar 

The use of the worldsheet system Scs, with its two types of coloured currents, v and v, 
is the reason for the Lagrange-multiplicr-type action (5.1). An even simpler example 
is the bi-adjoint scalar theory, BS in table 3. In this case, we can easily apply the 
procedure of [37] and obtain explicitly the equations of motion. As in that paper, which 
was concerned with the Einstein theory, the spacetime background fields modify the 
worldsheet theory only through the constraints. The deformation of the constraints in 
the bi-adjoint scalar theory is particularly simple: the deformed ambitwistor constraint 


becomes 


n = P 2 -> = P 2 + $ aa 'v a v ,a '+(t) aa 'v a v ,a ', (5.2) 


where we introduced currents for each of the two independent groups SU(Y) and 
SU(Y'). The equations of motion are obtained as anomalies obstructing the vanishing 
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of the constraint at the quantum level, 



+ simple pole. 


(5.3) 


If the equations of motion hold, there is no double pole and in fact the OPE is finite, 
because there can be no simple pole in the self-OPE of a bosonic operator in the absence 
of higher poles. The spacetime action associated to these equations of motion takes the 
Lagrange-multiplier form 



(5.4) 


It should be seen as the analogue of the gauge theory action (5.1). 


6 Discussion 

There are many issues to explore further and we briefly mention a few of them here. We 
have listed the central charges for the various models that are not already anomalous 
in table 4. It can be seen that many of the models have some critical dimension 
where the central charge vanishes. Indeed, one can often simply add some number 
of Maxwell fields to make them critical if one starts in low enough dimension. This 
suggests that a number of these models might give rise to plausible string expressions 
for corresponding loop integrands such as given in [12] for the type II theory in 10 
dimensions. However, an independent criterion is that the loop integrand so obtained 
should be modular invariant and this may well exclude many of the critical models as 
it does in conventional string theory. 

There is also the question as to whether there are further vertex operators that 
we have missed and therefore further sectors of these theories. For the 10 dimensional 
models, following [12], one can introduce a spin field 0“ associated to each T field and 
use these to introduce further vertex operators that will correspond to space-time fields 
with spinor indices. For the type II Einstein theory these give rise to the Ramond sector 
vertex operators [12] and it can be seen that the same procedure can be applied more 
generally to some of the models here, particularly the Einstein T*YM models. Following 
the same procedure one then extends the Einstein NS sector to include the Ramond 
sectors of type II gravity theories. However, we can see that the T*YM vertex operators 
can only be extended in this way on the one side corresponding to the spin operator 
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Theories 

Integrated vertex operators 

Central charge c 

E 

V h = (e • P + k • (i ■ P + k • ie • i^j 

3 (d- 10) 

EM 

v h , y 7 

Y 7 = (k • it ■ p) (e • P + k • ie • i^j 

3(d-10 + f) 

EMS 

V h , Vy, Vy, V S 

Vs = (k ■ i t ■ p) (jc ■ ip ■ tj 

3 (d 10+ m +™) 

BI 

V B i = (k • iik • ® 2 ) (e • P + k • ie • i^j 

\ (7 d - 38) 

Galileon 

V G = (k ■ iik ■ i 2 ) (k ■ 4q k • ^ 2 ) 

00 

1 

DBI 

Vbi, Vs bi 

v Sm = (k-iik-i 2 ) (k-it-pj 

5(7 d + m — 38) 

T* YM 

V a ={\t-\p,p]) ( e-p + k-m •*) 

Vg — (t {\p, p] + [q, y])) (e • P + k • 'he • 

|(c? — 12) 

ET*YM 

Vh, Vg , 

v g = (k-m-p+lt-\p,p])(i-p + k-m-&) 

Vg = (k-it- p + t- ([p, p] + [q, y})) (i-P + k-ie-i) 

3 (d- 10) 

NLSM 

v = (\t ■ [p,p]) (k ■ 'ilk • T 2 ) 

V = (t • ([p, p\ + [q, y})) (k ■ iik • § 2 ) 

3d - 19 


Table 4. Table of the different theories and their integrated vertex operators. 


constructed from the T in the Yang-Mills vertex operator. Thus one supersymmetry 
acts trivially on the Yang-Mills and hence is degenerate (it does not square to provide 
the Hamiltonian on the Yang-Mills fields). 




By extending the worldsheet matter fields we have generated new possible couplings 
to space-time fields. It would be interesting to explore whether these couplings can be 
made consistent in the fully nonlinear regime as described in [37, 38]. 

There remain other formulae based on the scattering equations, for which an un¬ 
derlying ambitwistor string theory has not yet been found. It would for example be 
interesting to find ambitwistor strings that give rise to the class of formulae with mas¬ 
sive legs [39-42], and that for ABJM theory [43, 44], although see the twistor string 
[45], 

Perhaps the most irritating issue is that we have not been able to find an Einstein- 
Yang-Mills model that is anomaly-free without unwanted linearized modes. Conven¬ 
tional string theory produces such amplitudes in open string theory and in closed string 
heterotic models. However, the ambitwistor heterotic string has corrupt gravity ampli¬ 
tudes and so far there has been no ambitwistor analogue of open strings. Nevertheless 
the T*YM model is likely to make sense and provide the correct amplitudes at 1-loop if 
modular, although the pure gauge sector does not have loop amplitudes beyond 1-loop. 
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A Correlators for Sym 

Here we give the proof of theorem 1. In the main text several versions of the present idea 
are realized. We will demonstrate and prove the mechanism in the simplest setting, 
which already contains all necessary ingredients, and comment on adaptations and 
restrictions afterwards. Concretely we use a single free fermion p a and a generic level 
zero current j a . The fields have the same OPEs as above, that is j a form a current 
algebra and p a are in the adjoint presentation of the j-algebra, i.e. 

P»P k (0) ~ V . j»/(0) ~ - , j»p‘(0) ~ i/“V ■ (A.l) 

a a a 

The strategy of the proof is as follows: both the full space-time amplitude A(g, h) 
and the world-sheet correlator A(g, h ) are a (multiple) sum of simple terms. The sum 
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in A is over trace sectors as well as a choice of gluon labels, while the sum in A(g, h) 
is simply the Wick expansion of the expectation value. Schematically we get 

A = A{x) and A = E A(y) (A.2) 

xzlX y€Y 

where A", Y are sets labelling the trace sectors and organization of sets of Wick con¬ 
tractions respectively. Then we will show that x E X x E Y and y E Y => y E X. 
Upon showing that each element in A", Y is unique we get X = Y. Along the way we 
will see that A{x) = A{x), hence establishing A = A. 

To clarify the structure of the discussion we firstly only insert integrated vertex 
operators on the world-sheet - which corresponds to considering the full Pfaffian in the 
CHY formula - keeping in mind that to get a non-vanishing result we need to go over 
to the reduced Pfaffian. That step will be taken at the end. 

So we will have to examine the correlation function of two types of operators, 

O gl = k-V t- p + t-j and O ar = A:-Te-T + e- P, (A.3) 


for (one half of) the gluon and graviton integrated vertex operators respectively. The 
claim is that the string-worldsheet correlator 


A(g,h) 



(A.4) 


where g and h are the sets containing the gluon and graviton labels respectively, is 
equal to (one part of the CHY representation of) the full space-time amplitude 

A= PfH, (A.5) 

trace sectors 


where the sum goes over all trace sectors possible. In particular, it includes a sum 
over the number of traces m — 1, • • ■ , [|g|/2]. The matrix n, defined in [22], of course 
depends on the trace sector. 

The main step in going between the representations two is the identity eq. (4.2), 
which we repeat here for the readers convenience 


Oab C(T) = fC(b,a\T) , 


(A.6) 


with /C, the ‘comb structure’ defined in the main text. Its arguments are the unordered 
set T and two of its elements, a,b E T. Using the anti-symmetry and multi-linearity of 
the Pfaffian, expression eq. (A.5) can be brought into the form 

E E !C(a u b 1 \T 1 )---lC(a m ,b m \T rn ) Pf M(h, {a,}, {b^h) . (A.7) 

trace ai<&i£Ti 
sectors 
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This is the representation of the amplitude which the world-sheet correlator eq. (A.4) 
will land us on. 

Let us now consider evaluating the correlator A(g, h). We will see that it gives rise 
to a multiple sum over terms, which turn out to be the same that eq. (A.7) sums over. 
The first step is to expand the product of all the O gl s into a sum. The sum is over all 
ways of putting either a kA> p or a j at each gluon insertion. This is a binary choice so 
it leads to 2^1 terms. Name the set of gluon labels which carry a k^p insertion by e 
for each term. The path integral over the T field can now be performed for each term 
individually. Since T is fermionic, the path integral vanishes unless \e\ is even. Define 
m := |e|/2, which is integer. The result of this path integral is of course simply a factor 
of 

Pf M(h, e\h) (A.8) 

for each term in the sum, by the standard reasoning described for example in [10] . note 
that, since e only appears once, the Pfaffian depends on the ordering of the elements 
in e. Now the correlator A is a sum over ways of partitioning g into e and g — e, with 
the condition that |e| be even, and each term in the sum looks like 6 

< id n j a ) PiM(h,e\h). (A.9) 

aEe aEg—e 

It should be clear that the remaining worldsheet correlator will give rise the product 
of /Cs and the remaining sum. Let us see how this happens in detail. Performing the 
Wick expansion of the p,j correlator breaks it down into a product of smaller pieces, 
so far until each factor contains precisely two (i.e. a pair of) p insertions accompanied 
by some subset of the j insertions. Label the pair of p insertions int the i th factor by 
a,i, bi and the accompanying set of j insertions by T % . Wick expansion makes sure that 
each choice of pairs and each choice of accompanying j insertions appears at least once 
and only once. Schematically we get 

m 

(x\pa n a* ) = sz p < Pai(>b ' n a* ) ■ ( A - iQ ) 

aEe aEg—e partitions i= 1 CiETi 

The remaining correlator is now easily evaluated using the OPEs eq. (A.l) to give 

( Papb JJ j c ) = lC(a, b\T) . (A.11) 

ceT 

c^a,b 

Note that the symmetry properties of the function /C in its arguments naturally arise 
from the statistics of the fields p,j. 

6 From now onwards we omit the colour structure and abbreviate t a ■ p(cr a ) = p a and t a ■ j(cr a ) = Jo- 
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Actually, performing the Wick expansion in eq. (A. 10) does not preserve the order 
of the p insertions, so, as they are fermions, a factor of (—1) might appear. We can 
absorb this factor by bringing the rows/columns of the matrix M into the same order 
as the p appear on the rhs of eq. (A.10). Then eq. (A.9) becomes 

^ K{ai,bi\Ti) ■ ■ ■ IC(a m , b m \T m ) Pf M(h, {aq, b x , ■ ■ - , a m , b m }\h) , (A.12) 


which is precisely the summand appearing in the full space-time amplitude 7 . We re¬ 
peat that Wick expansion ensures that every possible configuration of the summand is 
summed over, each term appearing at least once and only once. 

We have shown that the expressions A and A are sums over the same simple terms, 
involving /Cs and the corresponding Pf M. To clarify, on one hand, the sum in A goes 
over different ways of choosing m pairs a;, bi from g and different ways of forming m 
unordered sets T t from the labels left over, as well as the sum over m. The set X 
mentioned above contains as elements the ways of making such choices. On the other 
hand, the sum in A goes over ways of splitting the labels g into m unordered subsets 
Ti and the ways of picking a pair from each subset, as well as the sum over m. The set 
of these choices is Y. What remains to show is that the sums are actually the same 
or, equivalently, that each sum includes the other. To do so, we go back to the full 
expressions 


A(g,h)=i n n°2 r 


aeg 


aEh 


and A(g, h ) = C\ ■ ■ ■ C m Pf II . 

trace sectors 


(A. 13) 

and argue that if a term appears in A it also appears in A and vice versa. Additionally 
we argue that each term appears at least once and only once in each expression, which 
will conclude the proof that they are equal. 

It is clear that both sums contain the summation over m — 1, • • ■ , [|g|/2] in them, 
which is to be understood as the number of traces. Take a contribution from A with 
\e\ = 2m. Each term in the sum is uniquely determined my specifying m pairs {di,bi} 
and m unordered sets Ti. As mentioned previously, Wick expansion guarantees that 
each term appears once and only once. A given term should have a partner in A at m 
traces. Looking at the representation eq. (A.5) this is not straightforward to see, but 
the equivalent representation eq. (A.7) makes this readily apparent. The sum over trace 
sectors will include one term where the T* in A are precisely 8 the T) in A whereupon 
the sums a*, bi G T) will contain one term in which all a*, bi in A agree with those in 


7 In fact there will be additional sign factors from permutations the rows/columns in the Pfaffian. 

8 In a slight abuse of notation, what is called Ti in A is actually Ti U {cu, bi} in A. 
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A. This shows that each term in A has a partner in A , establishing the statement 
y E Y => y E X. Of course Wick expansion guarantees the uniqueness of the terms in 
Y. 

Conversely, one term in the summation in A is uniquely specified by fixing a trace 
structure and picking out one term of the summations over oq, b t . In other words, it is 
specified by a collection of m sets T* and a choice of pairs {a*, for each set. To see 
that any such term is also contained in A simply notice that the above data uniquely 
specifies a term in A via 

m 

IK***. II*> p f M {h,{{oi,bi}i}\h) . (A. 14) 

i =1 CiGTi 

Hence, each term in A has a partner in A and this establishes the statement x E X =>■ 
x E Y. The uniqueness of each element follows by construction. 

A.l The Reduced Pfaffian 

The Pfaffian we discussed so far actually vanishes for physical systems, i.e. when 
momentum conservation, gauge invariance and the scattering equations hold. Hence it 
is replaced by the reduced Pfaffian Pf n defined in either of the following equivalent 
ways 

( _No ( _W ( _ \a+b 

Pf'n := Pfn^v = A^ P fn a . = -> ; pm aj -, = Y^pfn a , b (A. 15 ) 

&a Oab 

where a, b label gravitons, with the restriction to not remove any row/column of the 
matrix B, and the i, j' label traces. On the other hand, we know that the expectation 
value of all integrated vertex operators will also vanish, and we have to insert precisely 
two fixed vertex operators. For an all graviton amplitude, this was discussed in [10]. It 
follows from BRST invariance that the amplitude is invariant under the choice of which 
vertex operators to take fixed/integrated. Hence, if there are at least two gravitons and 
arbitrarily many gluons, the full amplitude must be equal to the CHY formula. We will 
now show that the reduced Pfaffian also follows when using fixed vertex operators for 
two gluons or one gluon and one graviton, trying to present the following expressions 
in a suggestive form. 

A. 1.1 Two Gluons Fixed 

Denote the labels of the fixed gluon operators as c, d. With the reduced Pfaffian defined 
as 

pf , n = Pfn M , (A. 16) 
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there are two cases, f = i or j' ^ i 
from the Pfaffian and we can write 

■■■ Ci ■■■ Pfn M , 

with the gluons c, d being members of the trace T The factor fits into the inter¬ 
pretation of [10] as ghost field correlator. Note that there is no sum over choices of 
pairs in T), instead the comb K, appears with fixed start/end points, corresponding to 
the insertion of fixed vertex operators for the gluons c, d. 

In the second case ( j' ^ i), name the traces such that c G T\ and d e T 2 . Now each 
term in the expansion of the worldsheet correlator will look like (omitting all irrelevant 
factors) 


. In the first case the trace T t is totally removed 
= K{c,d\T$ Pfn M / , (A. 17) 


— V/C(c,a|Ti)/C(d,6|T 2 ) Pf (a, b, ■ ■ •) = C{T X ) C(T 2 ) V Pf (a, 6, - 

acd ^ y a <* 

b£T 2 b£T 2 


C(Ti)C(T 2 ) V — Pffa.Y'ffw &,•••) 

z —rr . / z —* 


aSTi 


@cd 


b£T 2 


C(Tx)C(T 2 ) V — Pf (a,-V^V 

aeT1 acd feeTi 


C(Ti)C(T 2 ) 


aGTi 

6 GTi 


&ac&db 

&cd 


Pf K &,•••) 


= C(Ti)C(T 2 ) y Pf (a, &,■■■) 

a<beTi Gcd 

= C(Ti)C(T 2 ) y (7 6a Pf(a,6,---)=C(Ti)C(T 2 )Pfn 2i2 /. 

a<bSTi 

(A. 18) 

Note that we had to use the scattering equations and the antisymmetry of the Pfaffian 
to arrive at the final result. 


A.1.2 One Gluon, One Graviton Fixed 

The computation for fixing one gluon and one graviton vertex operator is largely anal¬ 
ogous to the previous one. Moreover, BRST invariance guarantees that the final result 
will be as desired. Let us nevertheless demonstrate the necessary manipulations. De- 
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note the fixed gluon by c, with c £ T), and the fixed graviton by m 


— V/C(a,c|Ti) Pf(a, • • • ,m,---) = C(7\) V — Pf(a,--- ,m,---) 

amc aeTx aeT! ^ 

= C(T i)— Pf (V] CT ca a, ■ ■ ■ ,m, • ■ ■) 

amc ae3\ 

= C(Ti)— Pf (-<r cm m, ■ ■ ■ ,fh, ■■■) 

®mc 

= C(T, j Pf (m, - ■ , m, - ■) = C(Ti) Pfn, ,, . 


(A.19) 


Again we had to make use of the scattering equations. 


A.2 Adaption and Restriction 

As mentioned in the text, it seems not to be possible to find a level zero current via 
descent from p in a consistent way satisfying eq. (A. 11). Hence, the main text contains 
an adaption of the system discussed above, using two fermions p a ,p a , conjugate to each 
other. Via the descent, p a gives rise to j a while p a gives rise to j a . The OPEs between 
the currents and the fields are 


p“M/( o) ~ !/“V. 

z 

p“M?(o) ~ -/“V. 

z 


o) ~ , 

z 

p a (z)j c ( 0) ~ o . 


(A.20) 


We shall now examine the correlators of this system. 

First, note that by taking the fixed vertex operators to be (p + p), the discussion 
above would carry over verbatim. There is a crucial difference however. The current 
appearing in the associated integrated vertex operator does not quite satisfy eq. (A. 11), 
but instead 

( (Pa + Pa) (p b + p b ) n tic + Jc) ) = \T\ JC(a, b\T) . (A.21) 

ceT 

So each contribution from a different trace sector will come with a different prefactor 
nr mi, spoiling the relative coefficient between partial amplitudes. As the prefactor 
depends on the given partition of particles into traces, it cannot be removed by a field 
rescaling. The origin of this factor can be understood by simply counting the ways in 
which a full comb can be generated. If we represent the fields p, p by +, — and the 
currents j,j by ±, ++ respectively, the possible contractions can be found by drawing 
all allowed charge flows as in figure 2 

Observe that each contraction must have exactly one insertion of j (represented by 
++) or p (represented by +), independent of the length n of the chain, while there are 
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Figure 2. Charge flows 

n— 1 insertions of j or p. Summing over the possible positions of the tilded operator in 
the chain gives rise to the over-counting by \T\. Note that each contraction contributes 
exactly the same analytical & colour structure. 

Having understood the (non-trivial) origin of the factor \T\, the remainder of the 
discussion, showing how to remove it, follows trivially. Denote v the vertex operator 
containing p and j and v the one containing p and /, either integrated or fixed. It is 
now clear that choosing to insert v at m of the gluon punctures and v at the others will 
give rise (following the general discussion above) to the complete color ordered partial 
amplitude with m traces 

C{T 1 )---C{T m ) Pf'n , (A.22) 

which concludes the discussion. 
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